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ABSTRACT 

We study the Lax-Friedrichs scheme, approximating the scalar, genuinely 

nonlinear conservation law u t + f x (u) = 0, where f(u) is, say, strictly 
'* • 

convex, f > a* > 0. We show that the divided differences of the numerical 
solution at time t do not exceed 2(ta*) . This one-sided Lipschitz 

boundedness is in complete agreement with the corresponding estimate one has 
in the differential case; in particular, it is independent of the initial 
amplitude in sharp contrast to linear problems. It guarantees the entropy 
compactness of the scheme in this case, as well as providing a quantitive 
insight into the large-time behavior of the numerical computation. 
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Introduction 


We study monotonicity preserving schemes of the 3-point conservative form 
(1.1) v v (t+k) = v y (t) - A[h(v v (t),v v+1 (t)) - h(v v _ 1 (t),v v (t))], 

serving as consistent approximations to the scalar conservation law 
(1.2a) (x,t) +|| (u(x,t)) = 0 

and subject to the initial data 

(1.2b) v y (t) = u(x ,0), u(x,0) e D L DBV. 

t=0 

Here, v y (t) = v(x y ,t) is denoting the approximation value at the gridpoint 
(x v = vAx,t), k and Ax are, respectively, the time-step and mesh size such 
that the mesh ratio ^ = k/Ax is being kept fixed, and h(*,*) is the 
Lipschitz continuous numerical flux consistent with the differential 
one, h(v,v) = f(v). 

Studying conservative difference approximations to (1.2), one aims at 
having 

(i) compactness 
(ii) entropy condition. 

By compactness we merely mean the compactness of the family of solutions 

{v(»,t) = v(*,t;Ax), 0<t<T, 0<Ax<e}. A standard tool being used in that 

direction, e.g., [1], [3], [6], [11], is to guarantee that the total variation 

TV[v(t)] = I !v v+1 (t)-v v (t) i remains bounded in time, v e L (BV, [0,T]): 

since the mean value v(t) = ^ v y (t)Ax is independent of t, it then follows 

v 
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00/ 1 00 X 

that v £ L [L M L , [0,T]J; a classical argument which involves Helly s 

theorem, the diagonal process and Lipschitz continuity of |v(*,t)| ,, 

1 L 
implies L^^-compactness . By compactness there follows the existence of a 

subsequence Lj- limit solution v(x,t), v(x,t) = lim v y (t;Ax'), 

x=v>Ax' , Ax'+O 

0<t<T, satisfying (1.2) in the weak sense. In Section 2 we show that 3- 
point monotonicity preserving schemes are exactly those whose total variation 
does not increase in time; in particular, therefore, they admit a limit-weak 
solution. 

It is well known that independently of the initial smoothness, weak 
solutions of (1.2) are, in general, not unique. By the entropy condition, we 

refer to a variety of criteria which single out the so-called physically 

relevant solution, thus guaranteeing uniqueness: geometrically they require 

characteristics to propagate toward shocks; analytically they indicate the 
existence of vanishing viscosity. In the case f is convex, for example, 
they amount to Oleinik's (E) condition requiring u i e f t > u r ight across shock 
discontinuities. A standard way being used to verify the entropy condition, 
e.g., [1], [7], [13], is by constructing a discrete entropy pair satisfying an 
entropy inequality. Unfortunately, the limit solutions of monotonicity 
preserving schemes are not necessarily the physically relevant ones — 
examples of limit solutions violating the entropy condition in this case, are 
well known (e.g.. Example 2.4 below). 

Monotone schemes — those for which the RHS of (1.1) is non-decreasing in 

each of its v-arguments — is by now a classical example for a subclass of 

monotonicity preserving schemes, capturing both the compactness and the 
entropy requirements. This has been shown by successfully implementing ideas 
along the above lines, e.g. [1], [7], [13]. Recently, Osher [13] introduces, 
for the method of lines, a general E class of monotonicity preserving schemes 


\ 
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enjoying both properties of entropy satisfying compactness. 

Identifying 3-point conservative scheme according to their numerical 
viscosity coefficient [6], we arrive at Section 2, to the following concise 
characterizations! while monotonicity preserving schemes (compactness) are 
exactly those having numerical dissipation no more than Lax-Friedrichs (LF) 
scheme, no less than Courant-Isaacson-Rees scheme, entropy satisfying schemes 
are those further restricted by having no less dissipation than Godunov's 
scheme; the latter is, in fact, the fully discrete analogue of Osher's E 
s chemes . 

In both cases, the LF scheme plays a special role as the one having the 
most allowable numerical dissiptation. In Section 3 we begin studying the LF 

scheme in the genuinely nonlinear case where f is, say, strictly convex, 

. 1 

a* = Min f(v) > 0. We show that the following one-sided Lipschitz condition 
holds 


(L) 


sup 

V 


V v+I (t) “ v v-l (t) 
2 Ax 


ta. 


The one-sided Lipschitz bound is in complete agreement with the corresponding 
estimate one has in the differential case [10, Theorem 3.1]. In particular, 
it is Independent of the initial amplitude; this is in sharp contrast to the 
situation in the linear problem. Estimate (L) guarantees both compactness and 
the entropy condition as well as providing quantitive insight into the 
behavior of the numerical solution. In Section 4 we conclude, in verifying 
the (L) estimate for LF scheme when the CFL condition is replaced by the 
weaker monotonicity preserving requirement. 


*Here and elsewhere in the paper, * 
stands for f(v). 


denotes differentiation and a(v) 
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To motivate (L), one differentiate ( 1 . 2 ) to find that along an\ 


characteristic 


— = a(u(x,t)) we have (u^) = -aCu^) ; the latter 


equation is dominated by the characteristic ODE 


= -a*w , 


w = w(t) ■ 


For the last ODE one has w(t) < — - ; since physically relevant solutions of 

ta* 

(1.2) are exactly those whose characteristics can be drawn backward to the 
initial line t = 0, we conclude that u < —7- (weakly), which is the 

ta* 

differential analogue of having (L). Interestingly, the very same equation 
which rules out the existence of (long-time) strong solutions, e.g. [8], [ 9 ], 
can be used to show the existence of a physically relevant weak one. 
(Alternatively, the following simpler geometric argument prevails: the 

straight characteristics issued backward from (x^,t) and (x2»t) meet the 
initial line t = 0, at x l “ ta(u(x2»t)) and X2 - ta(u(x2»t)), 

respectively; the requirement for these characteristics not to intersect 
yields, after little rearrangement 


a(u(x 2 ,t)) - a(u(x, ,t) ) , 

£ — < - , 

X 2 - t 


* 2 > ^i* 


We close by saying that most likely the one-sided (L) condition holds for 
other schemes — those in the E class are, of course, natural candidates. If 
shown, (L) will provide, in the genuinely nonlinear case, a unified 
alternative to the standard total variation boundedness /entropy inequality 
approach, in showing the entropy satisfying compactness. 
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2. Three-Point Mbnotonicity Preserving Schemes 

We start by considering 3-point schemes in an increment form 

v v (t+k) = v v (t) + C^. 1/2 Av v - C~_ l/2 ^ v _ v 

C t+l/ 2 = c± (v v (t),v v+ i(t)). 

Having the scheme in such increment form is not a restriction. In fact we 
have 


(2.1a) 

where 

(2.1b) 


Lemma 2.1 

Every 3-point conservative scheme (1.1), can be written in an increment 
form (2.1) . Conversely, any 3-point increment scheme is conservative provided 
the following consistency requirement 


( 2 . 2 ) 

holds. 





v 




Proof : We follow Harten [5, Section II. 3] . Suppose (2.1) admits a 

consistent conservative form; equating the RHS of (2.1a) and (1.1) we obtain 
(for simplicity, we drop the time dependence) 

< + 1/ 2 Av v " C v-1/ 2 Av v- 1 = - X [ h(v v> v v+1 ) " 

Setting v v-1 = v y gives 
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(2.3a) 

while putting 
(2.3b) 


„+ _ , h < v v> v v> “ h(v v ,v v+1 ) 

C v+ l/ 2 “ X MT 


v v+1 , implies 


<f i. = X 

v- l/ 2 


h <v v ,v v ) - h(v v-1 ,v v ) 


Av 


v-1 


Thus, the incremental coefficients are necessarily those given in (2.3). We 
note that using (2.3), the consistency condition h(v,v) = f(v) amounts to 
having (2.2). 

Now, suppose (2.2) holds and define the consistent numerical flux 


< 2 * 4a v ) h(v v’ v v+l ) - f( V - T C v+1/ 2 Av v* 

Making use of (2.2) we find 

< 2 ' 4 Vi> h(v v-i' v v ) ■ - X - «V - x Vi/Vi' 

which puts the RHS of (2.1a) in the consistent conservative form (1.1), thus 
completing the proof. 


An essential property characterizing the scalar differential solution 
operator, which is highly desirable to be retained in the discrete framework 
as well, is preserving monotone profiles . We have 

Lemma 2.2 

Three-point monotonicity preserving schemes are exactly those whose total 
variation is non-increasing . They are characterized by the set of 
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Inequalities 


<2 - 5) < + v 2 >0 - c wv 2 >0 ' 1 -< + i / 2 - c ;;hV »• 

Proof : Differencing (2.1a) we obtain 

(2.6) Av v (t+k) = C+ +3/2 Av v+1 (t) + (1-C^.i^- C; +1/2 )Av v (t) + C^_ ^ Av^ (t) . 
Setting v v-1 (t) = v y (t) = v y+1 (t), i.e., A v v-1 (t) = A Vy (t) = 0, we find 

Av v (t+k) = C^. 3/2 A Vv+1 (t). 

By monotonicity preserving, the sign of Av y (t+k) must agree with that of 
Av v+ i(t), and hence, C y + 3/ 2 E C+ ^ v v+l ,v v+2^ should be non-negative for 
arbitrary v v+l ,v v+2* T ^ ie ot ^ er inequalities in (2.5) follow likewise by 
setting A v y _ x (t) = A v y+1 (t) = 0 and Av y (t) = A v v+1 (t) = 0. 

Next, we follow Harten [6] in showing that (2.4) implies the non-increase 
in total variation; summing (2.6), we obtain 

TV[v(t+k) ] i l |4v v (t+k)| < l |C^ 3/2 I |Av v+1 (t)| 

+ 1 Id - < + i /2 - < + v 2 >l l%d)l + l iv^l |4v u _ l( t)|. 

Re-indexing the first and third summations we find, on account of (2.4), that 
the RHS does not exceed 


) 
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l < + l / 2 l 4 V t >' + l d - Cl/r c wV l 4 v v (t)l 

+ l C wVo l 4 v v< t )l ■ l l 4 v v (t > I* 

V Z V 

Since, on the other hand, non-increasing total variation implies monotonicity 
preserving [6], there follows the equivalence between the two and their 
characterization by (2.4). 


The consistency requirement (2.2) shows that there is only one degree of 
freedom in setting up the recipe of 3-point conservative scheme; letting 


(2 ’ 7) Q v +1/ 2 H C v+V 2 + C v+V 2 ’ 

then by averaging (2.4a y ) and (2.4b y ) we have 

< 2 - 8 > »<v v v + i> - 1 ( f <V + £ <%h> - i Viii 1 ’.)' 

The scheme (1.1) then recast into the form 

(2.9) v v (t+k) = v v (t) - j [f(v v+1 (t)) - f(v v-1 (t))] +| [ A (Q v _l /2 A %_ 1 (t))], 

which reveals the role Q plays as the numerical viscosity coefficient [6]. 
Noting, that according to (2.2) 


( 2 . 10 ) 


_± _ V? + ^v+ V? 

v+ 1/2 2 


C v+V? C v+ V? B 1 r. x x 

? 2 ^ Q v+ 1/2 X Av v * 
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We conclude 


Corollary 2.3 

A 3-point conservative scheme (2.9) Is monotonicity preserving and total 
variation non-increasing. If and only If, Its numerical viscosity coefficient , 
Q v+ 1/ 2 » satisfies 


Af 


( 2 . 11 ) 


x tarl < Vi/< l ' 


We turn now to consider few examples. 

EXAMPLE 2.4 : Murman's scheme [12] is the nonlinear generalization of the 

upwind Courant-Isaacson-Rees (CIR) scheme 

(2.12a) v y (t+k) = v v (t) - A[ Svf i/ 2 (f v+1 -f v ) + (1 - s^i/ )(f v - f^)] 


where 


(2.12a 7 ) 


S v+V 2 


Af v 

Av. 


Af 
\ 

Av. 


< 0 


> 0 


1 Af v 

Noting that s \>f 1/ = y (* “ sign its incremental form (2.1a) amount 


to 


, , Af Af.. . Af Af 


(2.12b) 
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and it admits the consistent conservative form (1.1) with 

<2-i2c) mvW - 1 ( f v + £ v+r 4 %). 


Its numerical viscosity coefficient is given by 


(2. 12d) 


Af 


Q CIR,v+V 2 “ 


and hence it meets the monotonicity preserving requirement (2.11) provided the 

Af v 

CFL-like condition A sup |v — | <1 is fulfilled. If this is the case, then 

v Av v 

by Lemma 2.2 and the previously argued compactness, the scheme admits a limit 

solution v(x,t) = lim v y (t;Ax'). It is well known that this limit 

x=vAx',A x '-»-0 

solution nevertheless may turn out to be a physically irrelevant one as shown, 
for example, by choosing initial data v y (t=0) = sign(v-l^); with f(u) = u^, 
the scheme admits this initial discontinuity as an "expansive shock" steady 
state solution rather than dissolving it as a rarefaction, the reason being 
its lacking of (i.e., zero) dissipation in this case. 


EXAMPLE 2.5 : Godunov's scheme [4] is determined by the numerical flux 

(2* 13 a) h(v v* v v+l ) = f K (x \H-V 2 ,k) )’ 

where u^(x,t) is the solution of the Riemann problem (1.1) with initial data 

X * X v+V 2 
x " X v+V 2 



(2.13a') 


u R (x,t=0) = 
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Consider, for simplicity, the convex case, f > 0. Except for the sonic 



(2.8), A ^Qy+l^AVy equals the difference between the sum f v + f v+l and 
twice the numerical flux h(v ,v,. . ). Hence the numerical viscosity 

V V-f-i 

coefficient in this case is given by 


(2.13b) Q g , v-h l/ 2 


f v+l - 2f <l> + f v 



if a(v v ) < a(v) = 0 < a(v vfl ) 
otherwise 


We remark that in the general nonconvex case, one has due to Osher [13, 


Section 2], h^.v^) = ^ 


Min 
< v < v! 


max 


min /max 


'vf V 2 V ' V ^V 2 


[sign(v v+1 - v v ) f (v ) ] where 


v^ {j = Min/Max(v v ,v v+1 ). Hence, the numerical viscosity coefficient of 


Godunov's scheme boils down to A times . Max 

min , „ max 


fv+l + f v ‘ 2f(v) 


i.e., 


v v+l/ 2 < V < v v+V 2 


Av 


(2.13b') 

^G,v+ V 2 


i A Max 


f v+l + f v - 2f(v) 


Av. 




Af x 

Av. 


if a(v) vanishes, 
if a(v) * 0, 


min . . max 

V V+V9 < V < V ^V, 


min . , max 

v v+V 2 <t s V v+V 2 


Thus, Godunov deviates from CIR scheme exactly where the latter fails — it 
introduces nonzero dissipation in the sonic rarefaction case. 


EXAMPLE 2.6: 


The Engquist-Osher (EO) scheme [2] is having numerical flux 
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(2.14a) 


h(v v’ v v+l ) = I^ f v + f v+l " ^ |f(v)|dv), 

V,, 


and hence Its numerical dissipation coefficient in the convex case f > 0, is 
given by 


(2.14b) 


1 V+1 

%o,v+y 2 =x ir f lf(v)|dv=< 

*• V V 


fv+l - 2f( X) + f V 

f*g 

Af 

A I A — “I 

A v 


if a(v v )*a(v v+1 ) < 0 


otherwise 


Thus, it deviates from CIR scheme in the sonic case, introducing even more 
dissipation than Godunov does in the sonic shock case — in fact, it treats 
the latter case, a(v v ) > 0 > a(v^j_j) as a "compressive rarefaction. " 


EXAMPLE 2.7 ; The Lax-Friedrichs (LF) scheme is given by 

(2.15a) v v (t+k) = j (v v-1 (t) + v v+1 (t)) - j - f ( v v _ 1 (t)). 

Its incremental coefficients, numerical flux and numerical viscosity 
coefficient are given respectively by 

< 2 - 15t > < 4 v 2 -i T x ^- 

(2.15c) h<v V v vH-l ) “ 1 ( f v + £ vfl ' X 4v v)- 

Q LF,v+ l/ 2 s 


(2.15d) 


1 
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Thus, the LF scheme is monotonicity preserving provided the CFL-like condition 
Af v 

X sup |-t — | <1 is met, and hence admitting a limit solution in this case; 
v v v 

interestingly as we shall see later on, this solution still may be the 

physically irrelevant one. Only upon strengthening the CFL condition, 

X sup |a(v)| < 1, will we get the desired convergence due to the scheme 
v 

monotonicity, as is the case with Godunov and EO schemes. 

In view of the above examples, we see that the CIR and LF schemes have, 
respectively, the least and the most numerical dissipation allowed under the 
monotonicity preserving requirement, as (2.11), (2.12d), and (2.15d) read 


Corollary 2.8 

A 3-point scheme is monotonicity preserving if, and only if, its 
numerical viscosity coefficient , Q v+ i^ . satisfies 

^ 2 ‘ 16 ^ Q CIR,V+ V 2 < V 2 < ^LF , v+ V 2 * 

In [13], Osher introduces, for the method of lines, a class of E schemes 
which according to our terminology can be interpreted as exactly those having 
no less numerical dissipation than Godunov's; Osher showed that such E schemes 
satisfy the entropy condition. Carrying on his ideas to the fully discrete 
case we formulate though not prove 

Corollary 2.9 

A 3-point scheme is monotonicity preserving (compactness) and entropy 
satisfying if its numerical viscosity coefficient is further restricted by 
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(2.17) 


% , v+ 1/2 < V 2 < ^LF, v+ V 2 * 


Finally, we note that all the above mentioned schemes are first order 
accurate. This is not a coincidence, since generically we have 


Lemma 2.10 

Any 3-point monotonicity preserving scheme is of first order accuracy . 

Proof : The truncation error for smooth solution u of (1.2a), is given, 

modulo third order terms, by 

2 

- ~ j p — [ (Q(u,u) - X 2 a 2 (u)}u x J x . 

By Corollary 2.3, A|a(u)| < Q(u,u) < 1, which implies that the coefficient 
inside the inner curly brackets is non-negative 

Q(u,u) - A 2 a 2 (u) > A|a(u) [ • (l - A|a(u)|) > 0; 

Af v 

indeed, unless A |— — | = 1 where the scheme reduces to the trivial nongeneric 

v v 

case of pure translation, strict inequalities hold, excluding more than first 
order accuracy. We remark that an alternative proof can be given by the same 
argument used to show first order accuracy for monotone schemes: the first 

two inequality characterizing monotonicity preserving in (2.5) express 
according to (2.3), the fact that while setting v = v y ^ = v y = v v+ j» the RHS 
of (1.1), H(v v i> v v ,v v+l^’ is non decreasing in its first and third arguments; 
the third inequality in (2.5) implies that the partial derivatives w.r.t. to 



these arguments are non-negative; noting that the dependence on the second 
argument is not as essential, the result follows along the lines of [7]. 
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3. Time-Decay in the Genuinely Nonlinear LF Scheme 
The Lax-Friedrichs scheme 


(3.1) 


v v (t+k) = 


v v+l (t) + v v-l (t) 
2 




is essentially a staggered-type scheme. To simplify the notations we 
introduce the staggered differencing operator A = A(2Ax) and correspondingly 
D = A; D y (t) abbreviates Dv v (t), D(t) = sup D v (t). 

We are interested in the time decay of the numerical solution v y (t) in 
the genuinely nonlinear case, i.e., when a(v) * 0 for all v, say when 
f (v) is strictly convex 

a* = Min a(v) >0. 
v 


Theorem 3.1 

Consider the LF scheme (3.1) under the CFL condition 


(3.2) \ sup |a(v v (t)| < 1. 


Then, for arbitrary d, d > 1, we have 


S(t) ‘ Max t 7^7 r ' ’ TTdk S(t ‘ 0)1 - 

(t+dk)a* 


(3.3) 
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Proof : By Induction, starting with t = 0 where the second term inside the 

maximum brackets is reduced to D(t=0). The general step follows by showing 


(3.4) 


D(t) < 


dl 

t + dk 


M, 


M = Max[ — f- , AxD(t=0)]. 
AXa * 


Differencing (3.1) gives 


(3.5) 


Av v (t+k) 


^ v v+l (t) + *V-l (t) 
2 




By the strict convexity of f 


(3.6) 


~ * ~ .2 ~ ~ 
a(v )Av + — r- Av < Af < a(v „)Av 
' v' v 2 1 v 1 v v v+2 v 


* .2 
- |4 V 


The RHS inequality follows by second order Taylor expansion around v v+ 2> the 
LHS, around v . Inserted into (3.5), (3.6) yields 


Av v (t+k) < 


1 - *a(v v+1 ) 


2 Av v+1 + 2 


1 + ^(v^) 


Av v _i - |Av v+1 | 2 +|Av v _ 1 1 2 ) 


or, often division by 2Ax 


(3.7) 


1 - Xa(v 


D v (t+k) < 


v+1 


> ~ 1 + Aa ( v v +l ) ~ 

- Vl (t) + 2 Lo v-l (t) - 


ka J 




The CFL condition implies that the sum of the first two terms on the right 
does not exceed Max[D y+ j (t) ,D y ^(t)]; the sum of the two terms inside the 
last brackets on the right is not less than Max[D v+ ^ (t) ,D y ^(t)] 2 . Hence, 
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ka . 


(3.8) D v (t+k) < Max[D v+1 (t),D v _ 1 (t)] - _ Max[5 v+1 (t) ,5^ (t) ] ' 


We distinguish between two cases: 

t + dk 


(i) Assume M < 


ka*dX 


The quadratic z(l-Yz) is monotonic 


increasing as long as z < l/2y. This is the case with the RHS of (3.8) 

~ ~ ka* 

viewed as such quadratic in z = MaxfD^ (t) »D V _ 1 (t) ] with Y = -j- , since by 
assumption 


z - Ha45 v+1 (t),D v _ 1 (t)] < < 4- = 2? 


ka J 


ka . 


Thus increasing the RHS of (3.8) by replacing Max[D v+ ^ (t) ,D y ^(t)] with its 


dX 


assumed upper bound — ; — — M, we find 

t + dk 


~ j > ka*dX 

D v (t+k) < t + dk M U " 2(t + dk) M ) ‘ 


Since, by definition, M is greater than 


dXa a 


, the brackets on the right do 


k 

not exceed 1 - - + - — ^ , and the last inequality yields 


V t+k > * t + dk _ t + dk^ < t + k + dk M 


dX 


4* - 1 - 

(ii) M > — . The quadratic z(l - Yz) Has a maximal value 

ka*dX 

. Viewing the RHS of (3.8) as such a quadratic in 
~ ~ ka* 

z = Max[D .(t),D . (t) ] with Y = — ~ , it does not exceed a maximal value of 

VtI V— 1 Z 


(recall, d > 1) 
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D v (t+k) < 


1 


t + dk 


2ka. ka. (t+k+dk) 


dX 

t+k+dk 


t + dk 
ka^dX 


dX 

t + k +dk 


M. 


This completes the proof . 

We remark that (3.2) requires verification of the CFL condition 

X sup |a(v)| < 1, only at the gridvalues v = v (t). In the next section we 

V , Af v, 

discuss the situation when the weaker CFL-like condition X sup |-r — | < 1 due 

v % 

to monotonicity preserving holds. We pause here to discuss a few implications 

of Theorem 3.1, excluding the rarefaction-free trivial case where by (3.7) 

D(t=0) < 0 Implies D(t) < 0 later on. 

The one-sided Lipschitz estimate (3.3) involves a free parameter d, 

d > 1 which can be chosen so as to minimize the Lipschitz bound. 

2 

Choosing d = — , see Figure 3-1 , we conclude 

ka A D(t=0) 



Figure 3-1 
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Corollary 3.2 

Assume the CFL condition (3.2) holds . Then, the LF scheme (3.1) 
satisfies 


(3.9) 


sup 

V 


V v+1 (t) - Vy-l^) 
2Ax 


< 2D(t=0) < 

ta*D(t=0) + 2 



Proof : 
that d 
Implies 


We have to verify that the above choice of d Is admissible, i.e., 

° — — — > 1, or ? X. Indeed, CFL condition (3.2) 

ka*D(t=0) a*Av(t=0) 


• rsj 

Xa* sup Av (t=( 


To draw global estimates from the one-sided Lipschitz condition (3.9), we 
have — due to the staggered nature of the LF scheme — to distinguish between 
the even and odd numbered gridvalues 


<(t) * 


v 2v (t)» 


v v (t) H v 2v+l (t); 


under the CFL assumption of Corollary 3.2, (3.9) amounts to 


(3.10) 


AVy (t) 



Av”(t) < 


4Ax 



Given a grid function w = {Wy}, its increasing (similarily decreasing) 
variation, TV^tw] (similarily TV”[w]) is defined by 

TV*[w] = ± l Aw v . 

±AWy>0 
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Thus, we have 


(3.11a) 


TV [w] = TV + [w] + TV - [w] , 


while, when restricted to the Interval , 


(3.11b) TV [w] | - TV~ [w] i = X" - w( Xoo ) - w(x_ 00 ). 

— . oo * oo J * —oo* CD 1 =S X 

V — OO 


In particular, (3.11a), (3.11b) imply 


(3.11c) TVIwI j j - 2TV*M| ,*[»(,..) **(*.)]• 

1 _0O " CO 1 OO ' oo 

According to this terminology, (3.10) yields 

(3.12) TV + [v'(t)] < l < -L- SU pp[v(t)] 

A V ;(t)>0 ta* ta* 

and similarly for v"(t). 

Consider now the Cauchy problem (1.2). There are two cases: 

(i) The periodic problem . Let P be the period. Depending on whether 

it consists of an even or odd number of grid points, say even, we have by 

2P 

(3.1), that the increasing variation of v'(t) does not exceed — - — ; because 

ta * 

of its periodicity, (3.11c) imply that v'(t) has a total variation which 
equals twice its increasing one. Hence 

(3.13) TV [v' (t) ] < 7 — . 

t Min a(v'(t)) 

Thus, the total variation of v'(t) per period tends to zero as t tends to 
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infinity in an inverse linear rate, ~t In particular, since the mean 

value of v'(t), v'(t), is conserved in time, we conclude that for t large 
enough 

(3.14) sup|v'(t) - v'(t=0)| < - - - — • 

v ta(v'(t=0)) 

The last two estimates are — apart from the unessential factor of two — in 
complete agreement with the corresponding estimates one has in the 
differential case [10, Theorem 4.3]. 

(ii) The pure Cauchy problem . Suppose the initial data v(x=x v ,t=0) 
are supported in a finite interval of length L. Applying (3.11c) over the 
support of v'(t), we find that its total variation equals twice its 
increasing one; using (3.12) to bound the latter we conclude 


(3.15) 


TV[v' (t) ] 



supp [v ( t ) ] . 


In the differential case, the solution support at time t expands like t ^2 , 
[10, Section 4], and (3.15) implies a total variation decay < t - ^ . In the 
discrete case, the support of v(t) does not exceed L + 2X and hence 


(3.16) TV[v' (t) ] < - 4 “ + “T 1 , 

Xa * ta * 

showing the boundedness of the total variation. We remark that in both the 
periodic and pure Cauchy problems, similar estimates like those derived above 
hold for v"(t) as well, since the difference between its end values, 
according to the above v' -estimates, is uniformly bounded. Therefore, in 
either case, we end up with a bound on the total variation which shows that 



22 


the one-sided Lipschitz condition (3.9) guarantees the desired compactness, 
which in turn leads to the existence of a limit solution 

v(x,t) = 1 im v v (t;Ax'). It also implies the entropy condition: 

x=vAx^ , Ax'+O 

taking the limit Ax + 0 in (3.9), it follows that v i e ft > v right across 

discontinuities. Thus, v(x,t) = lim v y (t,Ax) is the unique, 

x=vAx, Ax+0 

physically relevant solution of (1.2). The desired convergence of LF scheme 
is, of course, well known. The point made here was to show that by virtue of 
the one-sided Lipschitz condition (3.9), one can deduce both the entropy 
satisfying convergence as well as quantitive insight into the large time 
behavior of the numerical solution. 


4. Time Decay in the Genuinely Nonlinear LF Scheme (cont'd.) 

In this section we continue our discussion on the time-decay of the LF 
scheme 


(4.1) 


v y (t+k) = 


v v+l (t) + v v-l (t) 




v+l 


- f vJ 


in the genuinely nonlinear case, say, strictly convex 

case, a* = Min a(v) > 0, restricted by the weaker CFL-like condition 
Af v V 

X sup |-r — | < 1, due to monotonicity preserving, 
v v v 

Theorem 4.1 

Let d, d > 1 be arbitrary and consider the LF scheme (4.1) under the 


CFL-like condition 
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(4.2) 


, , Af v, „ t + k(d-l) 

x 7 'xT' 4 

v 


Then (3.3) holds, l.e. , 


(4.3) 


D(t) < Max[ , —f— D (t=0) ] 

(t+dk)a . * + dlc 


Proof : We repeat the induction step in proving (3.4) from Theorem 3.1. 

Differencing (4.1) we find, as before. 


(4.4) 


Av v (t+k) = 


Av v+l (t) + Av^t) x 


2 ^ Af v+1 “ Af v-l) 


By the mean value theorem 


(4.5) Af y = a(v v )Av y , 


a(v ) = / a(v (6))d9, v (0) = v+ 0A V . 
6=0 


Inserted into (4.4), one reads 


1 + Aa(v v-1 ) „ 1 - Aa(v ) _ 

(4.6a) A v „(t+k) = ^ Av v-l (t) + ? Av^ (t) , 


'v+1' 


or, after division by 2Ax 


1 + Xa(v L ) „ 1 - Xa(v ) „ 

(4.6b) D v (t+k) = ^ L- Dv _ i(t) + - ~ D v+1 ( t) 


We distinguish between three cases: 

(i) The shock-rarefaction interaction case : Assume 

D (t) • D v+1 (t) < 0. If D v+1 (t) > 0 then D v _ 1 (t) < 0 and from (4.6b) 
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D v (t+k) < 


l - M Vv+1 ) 


D(t), while if D v _^(t) > 0, then 


1 + *a(v v-1 ) „ 

from (4.6b), D^(t+k) < ^ D(t). In either case, 


D v+ i(t) < 0 and 
(4.2) implies 


D v (t+k) < 


1 + X |a(v ) | 

=-^i — D(t> (l + 


t + k(d-l) i 
t + k(d+l) J 


dX = dX 
t + dk t+k + dk 


(ii) The shock case . D y _^(t), ^^^(t) < 0. Then by (4.6b) 

V t+k > 4 0 < t + i\ dk «• 

(iii) The rarefaction case: Assume D v _^(t), D v+ j(t) > 0. Our purpose 

is to reproduce the recursive inequality (3.8) which in turn led us to the 
desired bound (4.3). We rewrite with the help of (4.5) 

(4 * 7) &v +1 " * f v-l = a( Vfl ) 4 v v+l " a( %-l^ v v-l 

= a( %+l )(2:v v+l " ^v-l* + ( a( Vn> ~ a K,_i))Av v _i, 

where using the mean-value theorem once more we find 

1 

(4.8) a(v y+1 ) " a(v v _i ) = / [ a ( v v+i( 0 )) " a C v v _ 1 (®)) ]d9 

0=0 

1 1 

= / / a[nv (0) + (i-n)v v .(e)]dn[(i-0)Av . + eAv^Jae. 

e =0 n=o L 

Inserting (4.8) into (4.7) and (4.7) into (4.4) we obtain, after division by 
2 A x 
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(4.9a) D v (t+k) = 


1 + Aa(v j) „ 1 - Aa(v ) _ 

2 ~~ Vi<t) + 2 ~~ B v+1 (t) 


r 1 f 1 • ~2 

- k / / a [ • • •] dTi ( 1— 0 ) d6 (t) 

0=0 n=0 


1 1 


-k / / a[-]dtl 8d9D ,(t) D (t). 

0=0 n«=0 


Similarily, if we rewrite with the help of (4.5) 


Af v+1 - Af v-1 = a ^ v v-l^ Av v+l - Av v-l) + ( a(v v+l> - a(v v-l ) ) Av v+l 




and express the difference inside the second brackets on the right in terms of 
(4.8), substitution into (4.4) gives us after division by 2Ax 


(4.9b) D v (t+k) = 


1 + Xa ( v v_i> ~ 1 “ ~ 

2 Vl (t) + 2 D v+1^ 


r 1 f 1 • ~2 

- k / / a [ • • •] dn 0 d0 DV_. (t) 

0=0 n =0 

1 1 

- k / / a [ • • •] dn(l-0)d6 D .(t) D (t), 

0=0 n =0 1 1 


Since according to our assumption (t) > D V+ ^ (t) < 0, the fourth terms on 
the RHS of (4.9a) and (4.9b) are non-negative, while the sum of their first 
two terms does not exceed Max[D v _ 1 (t) ,D y+ ^ (t) ] ; we therefore find 


1 1 


D (t+k) < Max [5 . (t) ,D , (t) ] - k / / a[***]dn(l-0)d0 IT (t), 

1 0=0 n=o 


l l 


D (t+k) < Max [5 . (t) ,D . (t) ] - k / / a [ — ] dn 0 d0 (t) . 

v 1 0=0 n=o 



26 


Since a [ • • • ] > a*, each of the double integrals on the RHS of the last two 
inequalities is greater than V 2 > implying 

ka^ 2 

D v (t+k) < MaxtD^jOO.D^jft)] 2~ ® v -l.(t) 

lc3 

D v (t4k) < MaxtD^^t),^^)] - ~T ^(t). 

The last two inequalities amount to having (3.8) and the proof is completed as 
before. 

The one-sided Lipschitz bound appearing on the RHS of (4.3) is the 
maximum between the two terms 


2 

(t+dk)a^ 


dk 

t + dk 


D(t=0), 


each of which involves a free parameter d, d > 1* For d » 1, the first 
term is becoming negligible, dominated by the second D(t=0) — the initial 
strength of rarefaction. With d > 1, on the other hand, the first term 
dominates, provided no strong rarefactions are present. This is certainly the 
case after quite some time, when possible strong initial rarefaction were 
dissolved. Observe that as d decreases, we are forced to use a stricter CFL 
limitation 


A sup 
v 


Af. 


Av. 


„ t + k(d-l) 
t + k(d+l) ’ 


so that possible strong rarefaction will be dissipated. Thus, there are 
essentially two situations: 

(i) There exist no strong initial rarefactions, D(t=0) = 0(1) (in 
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particular, with smooth initial data). Then we can take d » 1 easing the 

CFL limitation and giving inverse linear time decay due to the first term 
-1 

~t • 

(ii) Strong initial rarefactions are present, D(t=0) = 0(h ^). Then we 

Af 

should choose d > 1 with CFL limitation X sup | 1 < — which over a 

v Av t k 
v 

time period T ~ X (i.e., after 0(h“^) time steps) will cause these strong 
rarefactions to dissipate. Afterwards, we are back in the first situation 

where no strong rarefactions are present, and again, there is a time decay 

-1 

~t 

It is clear that by using the time-uniform CFL-like requirement 
Af 

X sup |“ — | < 1 - 6, 6 > 0, we get the same inverse linear time decay; indeed, 
v Av 

v 

2 

taking d > - 1 will do for the above analysis. It shows that taking the 

full CFL limitation, 6 ~ 0, correspond to taking large d which may delay 
dissolving the rarefactions if initially strong — only after n ~ 5 ^ time 
steps will the strong rarefaction fully dissipate allowing the inverse linear 
time decay to dominate. 

Note, in particular, that as d -»■ °°, allowing the extreme CFL-like 
A f v 

requirement X sup |— — | = 1 to be used, is not enough to dissolve strong 
v Av 

v 

rarefactions if present, despite the scheme being monotonicity preserving in 

this case; the estimate (4.3) reduces to D(t) < D(t=0), and the following 

example bears out its sharpness. 

EXAMPLE 4.2 : Consider the LF scheme (4.1) applied to the Burger's equation, 

f (u) = u 2 , with initial data v y (t=0) = 6 and X = 1. The extreme CFL- 
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Af v 

like requirement X sup | 1 = 1 is trivially satisfied, yet the initial, 

v Av y 

physically irrelevant spike will travel one mesh to the right at the time, 
without being dissolved. The LF scheme amounts in this case, to a pure 

translation, lacking the dissipation to cause any decay. 

The above example does not contradict the convergence of the LF scheme to 
the physically relevant limit solution due to its monotonicity, as much as it 
shows, as was mentioned above, the importance of using the stronger CFL 
limitation A sup|a(v)| <1 — the correct one to guarantee monotonicity. 


v 
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